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The phase diagrams of low density Fermi-Fermi mixtures with equal or unequal masses and 
equal or unequal populations are described at zero and finite temperatures in the strong attraction 
limit. In this limit, the Fermi-Fermi mixture can be described by a weakly interacting Bose-Fermi 
mixture, where the bosons correspond to Feshbach molecules and the fermions correspond to ex- 
cess atoms. First, we discuss the three and four fermion scattering processes, and use the exact 
boson-fermion and boson-boson scattering lengths to generate the phase diagrams in terms of the 
underlying fermion-fermion scattering length. In three dimensions, in addition to the normal and 
uniform superfluid phases, we find two stable non-uniform states corresponding to (1) phase sepa- 
ration between pure unpaired (excess) and pure paired fermions (molecular bosons); and (2) phase 
separation between pure excess fermions and a mixture of excess fermions and molecular bosons. 
Lastly, we also discuss the effects of the trapping potential in the density profiles of condensed and 
non-condensed molecular bosons, and excess fermions at zero and finite temperatures, and discuss 
possible implications of our findings to experiments involving mixtures of ultracold fermions. 
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I. INTRODUCTION 

The lack of precise control over standard condensed 
matter systems has hindered the development of ex- 
periments that can probe systematically the effects of 
strong correlations. However, the large degree of con- 
trol in atomic systems, has made them powerful tools 
for studying many condensed matter phenomena, and 
in particular novel superfluid phases [l], H, H, 0, H, @. 
For instance, a current research frontier is the study of 
fermion mixtures with population imbalance @, H & fic| . 
Since the population of each component as well as their 
interaction strength are experimentally tunable, these 
knobs enabled the study of the BCS to BEC evolution in 
population imbalanced two-component fermion superflu- 
ids 0, H, H, E3| • In contrast with the crossover physics 
found in the population balanced case HHHIIl, these 
experiments have demonstrated the existence of phase 
transitions between normal and superfluid phases, as well 
as phase separation between superfluid (paired) and nor- 
mal (excess) fermions as a function of population imbal- 
ance [m, 

m he mm mils • 

Motivated by these recent experiments, there has been 
intense theoretical interest in understanding the phase 
dia gram of population imbalanced mixtures [2lL [22 . l23l . 

H, E3, m Ei|. So far, an accurate description 
of such mixtures is only available in the weak fermion 
attraction limit, and it is yet to be developed for inter- 
mediate fermion attraction around unitarity, or for the 
strong attraction limit. Some progress has been made 
in the strong attraction limit, where Fermi- Fermi mix- 
tures were described as a weakly interacting Bose-Fermi 
mixture [H, [26|, [3(| [M|, however, the effective boson- 
fermion and boson-boson scattering parameters were ob- 
tained only in the Born approximation. Strictly speak- 
ing, the Bose-Fermi description is valid only in the strong 



attraction limit, but may also provide semi-quantitative 
understanding of the phase diagram close to unitarity. 
Thus, the main goal of this manuscript is to analyze the 
boson-fermion and boson-boson scattering parameters 
beyond the Born approximation for arbitrary mass ratio 
of fermions, and use the effective Bose-Fermi mixture de- 
scription to generate improved phase diagrams and den- 
sity profiles of Fermi- Fermi mixtures with equal or un- 
equal masses in the strong attraction limit beyond the 
Born and mean- field approximations [26|, [3(| HH, [H, HH . 

The main results of this manuscript are as follows. 
First, we analyze three- and four- fermion scattering pro- 
cesses and obtain the exact boson-fermion and boson- 
boson scattering lengths as a function of mass anisotropy. 
Second, we use the exact boson-fermion and boson-boson 
scattering parameters to construct the phase diagram for 
Fermi- Fermi mixtures in the strong attraction limit. In 
addition to the normal (N) and uniform superfluid (U) 
phases, we find two different non-uniform phase sepa- 
rated (PS) states: (1) phase separation between pure 
unpaired (excess) and pure paired fermions (molecular 
bosons), and (2) phase separation between pure excess 
fermions and a mixture of excess fermions and molec- 
ular bosons, depending on the fermion-fermion scatter- 
ing parameter. The phase boundaries are very sensitive 
to the masses of the fermions, and also to the boson- 
fermion and boson-boson interactions. For equal mass 
mixtures, our results for the phase boundary between 
the PS(2) and the uniform (U) phase improves on previ- 
ous saddle-point (mean-field) results, and the quantita- 
tive changes are substantial, but not dramatic. However, 
there is a dramatic increase in quantitative differences 
between mean-field and the present results for unequal 
mass mixtures as the mass ratio deviates from one. In 
particular, these differences are more pronounced when 
heavier fermions are in excess indicating the importance 
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of taking into account scattering processes beyond the 
Born approximation. Furthermore, we discuss the effects 
of the trapping potential on the density profiles of con- 
densed and non-condensed molecular bosons, as well as 
excess fermions at zero and finite temperatures. Lastly, 
we discuss the implications of our findings to possible ex- 
periments involving Fermi-Fermi mixtures with equal or 
unequal masses and equal or unequal populations. 

The remainder of this manuscript is organized as fol- 
lows. In Sec. HIl we discuss briefly the Hamiltonian for 
Fermi-Fermi mixtures with equal or unequal masses and 
emphasize that the system reduces effectively to a Bosc- 
Fermi mixture of molecular bosons and excess fermions 
in the strong attraction limit. In Sec. IIII1 we analyze the 
exact boson-fermion and boson-boson scattering lengths 
as a function of mass anisotropy, which are used to calcu- 
late the resulting phase diagram of Fermi- Fermi mixtures 
in the strong attraction limit. In Sec. IIV1 we discuss the 
stability of the effective Bose-Fermi mixture of molec- 
ular bosons and excess fermions in three, two and one 
dimensions. We analyze the stability of population im- 
balanced Fermi-Fermi mixtures in the strong attraction 
limit in Sec. IIVB1 and we construct the phase diagrams 
of these systems in Sec. I IV CI In Sec. [V] we discuss the 
effects of harmonic traps on the density profiles of con- 
densed and non-condensed molecular bosons, and excess 
fermions at zero and finite temperatures, and show their 
experimental signatures. Lastly, we present a summary 
of our conclusions in Sec. IVII 



II. HAMILTONIAN 

In order to calculate the correct phase diagrams of 
Fermi-Fermi mixtures in the strong attraction limit, it is 
necessary to obtain first the correct scattering parameters 
between two Bose molecules (paired fermions), and also 
between a Bose molecule and an unpaired fermion. To 
achieve this task, we begin by describing the Hamiltonian 
density for a mixture of fermions (in units of h = kg = 1) 
as 



H(x) = 5^(3 



2m„ 



V'o 



9^{x)^i{x)ijji{x)ijj^{x) 



(1) 



where ty a [x) is the field corresponding to the creation of 
a fermion with pseudospin index a, at position and time 
x = (r, r). Here, g > is the strength of the attractive 
fcrmion-fcrmion interaction and a identifies two types of 
(t and |) fermions. This notation allows the analysis of 
a mixture of fermions with equal or unequal masses, as 
well as equal or unequal chemical potentials. Through- 
out the manuscript, we assume that the lighter fermions 
are always f and that the heavier fermions are always |, 
as intuitively suggested by the direction of the arrows, 
and that the chemical potentials fix the population 
(density) n a of each type of fermion independently. 



The contact interaction Hamiltonian given in Eq. ([T]) 
generalizes the equal mass and equal chemical poten- 
tial Hamiltonian that is used to study the BCS to BEC 
evolution within the functional inte gra l formalism [l3j |. 
The functional integral formulation [13l l26l . [30L l3l| cap- 
tures some essential features of the evolution from BCS 
to BEC superfluidity for Fermi mixtures with equal or 
unequal masses as well as with equal or unequal pop- 
ulations. However, truly quantitative results are cur- 
rently possible only in the BCS limit, where the theory is 
simple, but the temperatures required to reach the BCS 
regime are very low and hard to be achieved experimen- 
tally. In the unitarity regime, experiments can be per- 
formed and the phase diagram can be explored since the 
critical temperature for superfluidity is attainable, but 
an accurate theoretical description of this regime is still 
lacking. While in the BEC limit, not only the tempera- 
ture required to reach the BEC regime is experimentally 
achievable, but also the theory becomes simple since the 
Fermi-Fermi mixtures can be described effectively by a 
weakly interacting mixture of molecular bosons and ex- 
cess fermions [22], [26|, [3(J HH . However, the initial propos- 
als of such effective Bose-Fermi mixtures [22, 26] can pro- 
vide only semi-quantitative results for comparison with 
experiments in the BEC regime since the scattering pa- 
rameters for two molecular bosons and a molecular boson 
and an excess fermion are obtained only in the Born ap- 
proximation. In order to overcome this defficiency, we 
discuss next the boson-fermion and boson-boson scatter- 
ing parameters beyond the Born approximation for arbi- 
trary mass ratio. The correct scattering parameters will 
be used in Sec. IIV CI to construct phase diagrams and 
density profiles for quantitative comparisons with exper- 
iments in the BEC regime. 



III. BOSON-FERMION AND BOSON-BOSON 
SCATTERING LENGTHS 

Mixtures of two types of fermions in the strong attrac- 
tion limit can be described by effective Bose-Fermi mod- 
els [22], [2^, [3(| [M| , where fermion pairs behave as molec- 
ular bosons and interact weakly with each other and with 
excess unpaired fermions. Scattering lengths between 
two molecular bosons (ass), and between a molecular 
boson and an excess fermion (bbf) were calculated in 
the Born approximation using many body techniques for 
equal [HI, [22j and unequal masses [26|, |3(|. However, 
these results do not agree with calculations using few 
body techniques [H, [3g|, because it is necessary to go 
beyond the Born approximation. 

In the case of fermion mixtures with equal masses, 
while the Born approximation in many body theory leads 
to clbb — 2a,p [l]| and clbf — 8ap/3 [H, [26J, the results 
from few body techniques are asB ~ 0.60a_p [H[ and 
asF ~ 1.18of [3^|. However, a diag rammatic approach 
beyond the Born approximation [38l 39] for equal mass 
fermions recovers the few body results. In this section, 
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we generalize these diagrammatic approaches and ana- 
lyze the boson-fcrmion and boson-boson scattering pa- 
rameters for two types of fermions with unequal masses 
in order to make quantitative predictions for experiments 
involving Fermi-Fermi mixtures in the strong attraction 
limit. Here, we show that the diagrammatic many body 
approach for unequal mass fermions produce results con- 
sistent with three and four body techniques that were 
recently used to obtain clbf and asB as a function of 
mass ratio [3^, [36] . We analyze the technical aspects of 
the boson-fermion scattering parameter for unequal mass 
fermions in some detail, since they are much easier to 
present, while we do not discuss in great detail the tech- 
nical aspects of the boson-boson scattering parameter for 
unequal masses, as they are extremely cumbersome. De- 
tailed descriptions of the boson-fermion and boson-boson 
scattering parameters for equal mass fermions can be 
found in the literature 38, 391. 
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FIG. 1: Diagrammatic representation of the integral equation 
for the boson-fermion scattering T-matrix Tjf F (p,po). Here, 
we use the notation (— |) =]. and vice versa. 

We begin our analysis by describing a zero tempera- 
ture (T = 0) diagrammatic representation for the boson- 
fcrmion scattering T-matrix Tjf (p,Po) as shown in 
Fig. [TJ where wf — k 2 /(2m F ) and wb — k 2 /{2ms) are 
the kinetic energies for the excess fermions and molecu- 
lar bosons, respectively, and = — l/(mna F ) < is the 
binding energy of the molecular bosons. Here, 



mg = TOj + mi 



m-f + mi 



(2) 
(3) 



are masses of the molecular bosons and twice the reduced 
mass of the ] and { fermions, respectively. In this figure, 
single lines represent retarded free fermion propagators 



G , (T (k, w) 



1 



(4) 



where w a = k 2 /(2m a ) is the energy and \x a is the chem- 
ical potential of the cr-type fermions. Similarly double 
lines represent the retarded molecular boson propagators 



D (k,w) 



A I 3 / 2 

47r/TO| , 
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(5) 



obtained from the expansion of the effective ac- 
tion (26l . [30j . and a corresponding RPA re-summation 
of the fermion polarization bubbles leading to 
-Do(k, w) = — g/[l + gT(k,w)] where the fermion 
polarization bubble is r(k, w) = J2q qo Go,a(k + 
q, w + (7o)Go,cr(— q, — go)- Integration over the internal 
momentum q and frequency go leads to T(k, w) = 
r(0, 0) + [m^ 2 /(47r)] [ Wg - w - ^ - ^ - i0+) 1/2 which 
in combination with the definition of the fermion- fermion 
scattering length clf = mnT FF (0, 0)/47r, and the 
fermion-fermion T-matrix T FF (0, 0) = -g/[l + gT(0, 0)] 
lead to the final result described in Eq. ([5]). 

On the right hand side of Fig. [TJ the first diagram 
represents a fermion exchange process, and all other (in- 
finitely many) possible processes are included in the sec- 
ond diagram. In all diagrams, we choose f Q) to la- 
bel lighter (heavier) fermions such that lighter (heav- 
ier) fermions are in excess when F =1 (F =|). This 
choice spans all possible mass ratios. In the following, 
we set [i a = since all of the calculations are performed 
for three fermions scattering in vacuum. The T-matrix 
(p,Po) satisfies the following integral equation 



rpBF 



T k (p,p ) = -G - F (k + p, w B - w F + e b + po) 

- ^2 A)(q, w B + £6 + ?o)G f ,ir(-q, Wf - qo) x (6) 
q,9o 

F (q, q )Ga- F {p + q, w B - w F + £fc + Po + qo), 

where we used (— f) =J. and vice versa. On the right 
hand side, we can sum over frequency go by closing 
the integration contour in the upper half-plane, where 
Tjf (q, qo) and Do(q,WB + e& + go) a re analytic func- 
tions of go- Since this integration sets go = k 2 /(2mp) — 
q 2 /(2mp), we set p = k 2 /(2mp) — p 2 / (2mp) in order to 
have the same frequency dependence for the T-matrix on 
both sides [3!| . Since we are interested in the zero-range 
low energy s-wave scattering, we average out the angular 
dependences of k and p. When k — ► 0, the generalized 
integral equation defined in Eq. ([7]) can be expressed in 
terms of the boson-fermion scattering function ak^o(p) 
as 

m U a o F (p)/ m BF 



1 



dp 1 + (mnp 2 /mBF + dp 2 ) 1 / 2 p 2 + ■ 

dq , ( q 2 + 2m F qp/m B + P' 
— In ' 
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BF 
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Here, we used the definition of the boson-fermion scat- 
tering length 



BF 



(P) 



m B F 
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mBF 



1/2' 



Trip), 

(8) 

with its full momentum dependence, where rngj? is twice 
the reduced mass of an excess fermion and a molecular 
boson given by 

2msmp 



rriBF 



mB + m,F 



(9) 
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The integral equation shown in Eq. {7J) as well as the scat- 
tering length expression shown in Eq. ([8]) reduce to the 
results for the equal masses |38l.l39l| when mj = m| = m. 
Notice that only the fermion exchange process is taken 
into account in the Born approximation, and that ne- 
glecting the second term on the right hand side of Eq. ([7]) 
leads to a BF (0) = 2(mBF/mn)aF which is consistent 
with our previous results [2^, [30]. However, we need 
to include both terms in order to find the exact boson- 
fermion scattering length. 
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FIG. 2: The boson-fermion scattering length a BF (0) / 'aj? ver- 
sus mass anisotropy m|/mj is shown in (a) when lighter j- 
type (hollow circles) or heavier (solid circles) j-type fermions 
are in excess. The boson-boson scattering length a BB (0)/aF 
versus mass anisotropy mf/m^ is shown in (b). In addition, 
the scattering lengths in the Born approximation are shown 
as solid lines. 

Next, we solve numerically the integral equation given 
in Eq. ([7]), and obtain a BF (p) as a function of the mass 
anisotropy m-p/m^. The exact solutions and the Born 
approximation values of a BF (0) are shown in Fig. HJa). 
When mf = mj_, we find a BF (Q) « 1.179a^, which 
is consistent with results previously found for equal 
masses [U [H H, HE HI 111 • Notice that a B F {0) de- 
creases (increases) from this value with increasing mass 
anisotropy when the lighter (heavier) fermions are in ex- 
cess. In addition, in the limit of m| /m^ — > 0, while 
a BF (0) — > aj? when the lighter fermions are in excess, 
a BF (0) grows rapidly when the heavier fermions are in 
excess. Notice also that the Born approximation values 
for a BF (0) are not in agreement with the exact values for 
any mass anisotropy, but the general qualitative trends 
are captured by the Born approximation as can be seen 
from Fig. Ufa). 



In addition, we present results for a BB (0) as a function 
of mass anisotropy i«|/m| in Fig. [2jb) . The exact results 
of the boson-boson scattering parameters for unequal 
mass fermions can be obtained either by extending the 
diagrammatic approach for equal mass fermions (3a . |39| 
or by using few body techniques [35| . The Born ap- 
proximation values a BB (0) = (tub /fn-\{)<iF are found 
in Ref. [26|, HI ■ Since the technical details to calculate 
the boson-boson scattering parameters are quite cum- 
bersome, and not particularly illuminating, here we just 
mention that the results found in the literature (35j can 
also be obtained diagramatically for the unequal mass 
case. As shown in Fig. [Hb), a BB (0) grows very slowly 
as the mass ratio m^jm^ decreases, in constrast with 
the much more rapid growth of the Born approximation 
values of a BB (0). As expected, the Born approximation 
values of a BB (0) are not in agreement with the exact 
values for any mass anisotropy. 
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TABLE I: Exact boson-boson (aBB) and boson-fermion (obf) 
scattering lengths for a list of two-species Fermi-Fermi mix- 
tures. Here, asi (ibi) corresponds to excess-type of | (1) 
fermions. 

Several atomic gases of fermions ( 6 Li, 40 K, 87 Sr [ipj . 
and 171 Yb [4l|) are being currently investigated, and ex- 
perimental methods for studying mixtures of two types 
of fermions are being developed in several groups. Thus, 
anticipating future experiments involving mixtures of two 
types of fermions, we show in Table [I] the boson-fermion 
and boson-boson scattering lengths for a few mixtures. 

Here, we would like to make two comments. First, 
it is quite remarkable that the diagramatic many body 
approach recovers the few body results for boson-fermion 
and boson-boson scattering lengths for arbitrary mass ra- 
tio m| /mi. The diagramatic many body approach is per- 
formed in momentum space, while the few body approach 
is performed in real space. The two methods are equiv- 
alent because all of the possible scattering processes are 
taken into account exactly in the diagramatic many body 
approach at T — for three or four fermions. However, 
the calculation of the scattering parameters for three or 
four fermions in the presence of many others (arbitrary 
number of particles) at zero or finite temperatures re- 
quires a full many body approach, which is inevitably 
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approximate and more difficult to implement. Second, 
our analysis does not include the effects related to Efimov 
(three body bound) states, which may become important 
when analyzing the scattering parameters as a function 
of mass ratio m-f/mi (36| . In particular, the mixtures 
of 6 Li and 87 Sr or 6 Li and 171 Yb have mass ratios of 
rriLi/msr ~ 0.068 and mLi/myb ~ 0.035, which are be- 
low the critical ratio m^/m^ rs 0.073 for the emergence 
of Efimov states. 

Having presented the boson-fermion and boson-boson 
scattering lengths for arbitrary mass ratio m-\ / mi , we 
discuss next the resulting phase diagrams for Fermi- 
Fermi mixtures in the strong attraction limit, where the 
system can be effectively described by a Bose-Fermi mix- 
ture [22|, HE [3(| HH of molecular bosons and excess un- 
paired fermions. 



IV. BOSE-FERMI MIXTURES AT ZERO 
TEMPERATURE 

In this section we use the effective Bose-Fermi mix- 
ture description [22|, Hfl [3l| to analyze the phase dia- 
gram of population imbalanced Fermi-Fermi mixtures in 
the strong attraction limit. We describe first the gen- 
eral stability conditions for Bose-Fermi mixtures at zero 
temperature, and use this connection to discuss the sta- 
bility and phase diagrams of Fermi- Fermi mixtures in the 
strong attraction limit. 



A. Bose-Fermi Mixtures 

The ground state of Bose-Fermi mixtures can be de- 
scribed by the free energy (3(J 52 



£ — £b j t£f 



Ubbu] 



~U BF n F n B LlF n F~ ^B n B > (10) 



which characterizes the center-of-mass degrees of free- 
dom for a mixture of single-hyperfine-state bosons and 
fermions. Here \xp and np (/is and Ub) are the den- 
sity and chemical potential of fermions (bosons), ep is 
the Fermi energy of the fermions, and Ubb and Ubf are 
the repulsive boson-boson and boson-fermion interaction 
strengths. The density of single-hyperfine-state fermions 
in three dimensions is given by np — (1/V) X^k F 1 = 
fc|n/(67r 2 ), where kp is the Fermi momentum and V is the 
volume. The first term in Eq. (TIT)]) is the total kinetic en- 
ergy of bosons, which is assumed to be much smaller than 
all other energies, and is neglected. This assumption is 
very good since essentially all bosons are condensed in the 
k = state, when the boson-boson and boson-fermion in- 
teractions are weak. The second term in Eq. (fit))) is the 
total kinetic energy of fermions, which in three dimen- 
sions is given by £p = (1/V) 5Zj| c ' <fcF ek ,- F = 3ei?n^/5, 
where £k,F = \k\ 2 /(2mp). 



From the free energy given in Eq. (|10[) . we obtain the 
fcrmion and boson chemical potentials using the condi- 
tion d£/dn.i = with i — F, B, leading to 



up = e F + Ubfub, 

Ub = Ubbit-b + Ubftif. 



(11) 
(12) 



Then, we use the positive definiteness of the Bose-Fermi 
compressibility matrix Kij = dfj,i/drij, 



du F du B _ da F dfi B 
dnp driB dns dnp 



>0, 



(13) 



and find that bosons and fermions phase separate when 
the condition 



4tt 4 U% r , s 

(3D) 



(14) 



is satisfied in three-dimensional systems [30J, |42J . There- 
fore, the stability of uniform superfluidity puts an upper 
limit on the density of fermions in three-dimensions. 

Following a similar approach in lower dimensions, 
where np = k F /(4:ir) and £p = epnp/2 in two dimen- 
sions, and np — kpjit and £p — epnp/S in one dimen- 
sion, we find that the bosons and the fermions phase 
separate when the conditions 



mp U 2 BF 
' j bb' 



~ 2tt U 



np < 



mp U BF 
Ubb' 



7T 



(2D) 
(ID) 



(15) 
(16) 



are satisfied, respectively, for two- and one-dimensional 
systems. Notice that, the stability of uniform superfluid- 
ity puts a lower limit in one dimension, which is in sharp 
contrast with the three-dimensional result. Furthermore, 
the stability condition in two dimensions does not depend 
explicitly on the density of fermions (see also Ref. [43j|). 
However, the results in lower dimensions have to be used 
with caution, since quantum fluctuations are more pro- 
nounced, and may affect these stability conditions. 

For an atomic Bose-Fermi mixture, we can also de- 
scribe analytically a finer structure of phases. There 
are four possible phases [42j: (I) PS(1) where there is 
phase separation between pure fermions and pure bosons; 
(II) PS (2) where there is phase separation between pure 
fermions, and a mixture of fermions and bosons; (III) 
PS (3) where there is phase separation between pure 
bosons, and a mixture of fermions and bosons; and (IV) 
PS (4) where there is phase separation between two dif- 
ferent mixtures of fermions and bosons. 

For a three-dimensional weakly interacting Bose-Fermi 
mixture, we follow Ref. [42j and find that there are only 
two stable phases within the phase separation region: 
(I) PS(1) where there is phase separation between pure 
fermions and pure bosons, and (II) PS (2) where there 
is phase separation between pure fermions, and a mix- 
ture of fermions and bosons. We obtain analytically the 
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condition 



n p > 



1125tt 4 U% B 
128™* ire F 



5 Ubb 
AU B f 



n B , (3D) (17) 



for the transition from the PS(2) to the PS(1) phase [301 ] . 

In lower dimensions, we find that the structure of the 
phase diagram is quite different. In two dimensions, 
the phase separated region consists only of PS(1) where 
there is phase separation between pure fermions and pure 
bosons. While, for a one-dimensional weakly interact- 
ing Bose-Fermi mixture, the phase separated region con- 
sists also of two regions: (I) PS(1) where there is phase 
separation between pure fermions and pure bosons, and 
(III) PS(3) where there is phase separation between pure 
bosons, and a mixture of fermions and bosons. We obtain 
analytically the condition 



riF < 



3m F U BF Ubb 
2tt 2 U bb ~ Ub~f~ 



n B , (ID) 



(18) 



for the transition from the PS(3) to the PS(1) phase. 
Notice that the structure of the PS (3) phase in one- 
dimension is very different from the structure of the 
PS (2) phase in three dimensions. Again, the results in 
lower dimensions have to be used with caution, since the 
quantum fluctuations are more pronounced, and may af- 
fect these stability conditions. 

Next, we concentrate only on the three-dimensional 
case, and use the stability conditions found above as well 
as the interaction (scattering) parameters obtained in 
Sec. IIIII to analyze the phase diagrams of Fermi-Fermi 
mixtures in the strong attraction limit. 



B. Fermi- Fermi Mixtures 
in the Strong Attraction Limit 

To make an analogy between Bose-Fermi mixtures 
and population imbalanced Fermi-Fermi mixtures in the 
strong attraction limit, we identify F = {f or J,} as the 
excess fermions. This identification leads to the density 
of excess fermions (tie) and molecular bosons (n-s) given 
by 



r>E 



Up — Tl-F 



HU 



n 



n B = 



n E 



n-F, 



(19) 
(20) 



respectively, where n = n-\ + rij is the total density of 
f- and |-type fermions. Here, we use (— j) =\ and 
vice versa. For instance, if F =| fermions arc in excess, 
the density of excess fermions and molecular bosons are 
he = Tif — ri[ and ns = (n — n E )/2 = ni, respectively, 
such that all |-type fermions are paired with some of 
the f-type fermions to form molecular bosons, but there 
are f-type fermions left unpaired. It is important to em- 
phasize that the internal degrees of freedom (electronic, 
vibrational, and rotational) of molecular bosons are not 



explicitly considered here, in the same spirit of the de- 
scription of atomic bosons presented in Sec. IIV Ai where 
the electronic degrees of freedom were also not explicitly 
considered. 

For three dimensions, we define the boson-boson and 
boson-fermion interaction strengths 



Ubb 
Ubf 



Ana b b 47T7i 
m B 

AirClBF 

m B F 



m B 

47T/3 F 

1 

m B F 



-of, 



-of, 



(21) 
(22) 



where a F , o-bb — Ib&f and olbf = Pfo-f are the 
fcrmion-fcrmion, boson-boson and boson-fermion scat- 
tering lengths. Here, 73 = ubb/clf and 0p — o,bf/<1f 
are constants, which are found in Sec. IIIII as shown in 
Fig. Hand Table |U In addition, we define the population 
imbalance parameter 



P = 



(23) 



such that \P\ — ns/n, and n = K F /(3ir 2 ), where N a 
is the number of cr-type fermions and Kp is the Fermi 
momentum corresponding to the total density of fermions 
defined by K% = (k% d + k F ] )/2. 

Using these definitions, the phase separation condition 
Eq. ([M]) becomes 



\P\ > 



^lB m BF 



(24) 



where A = 1/(KfCLf) is the scattering parameter. Simi- 
larly, the condition given in Eq. (|17|) becomes 



\P\ 



^ _ 5j B m B F \ > 33757r 3 7|m| f 3 
8/3 F m B J ~ 81920%m B m% 



5j B mBF 



^PfTUb 

(25) 

for the transition from the PS(2) to the PS(1) phase. 

We would like to emphasize that these analytic expres- 
sions given in Eqs. and (f2l)]) are valid only in the 
strong attraction limit when 1/(KfCLf) 1, but they 
may still give semi-quantitative results for 1/ '(K F a F ) > 
1. Close to the unitarity, the Bose-Fermi description of 
Fermi-Fermi mixtures in terms of molecular bosons and 
excess fermions is not reliable, since the binding energy 
of molecular bosons is small and the interactions between 
molecular bosons and excess fermions or between two 
molecular bosons may be sufficient to cause dissociation 
of the molecules into directly scattering fermions. How- 
ever, there may be an intermediate regime between uni- 
tarity and the strict BEC limit where we can describe 
Fermi-Fermi mixtures in terms of a mixture of molec- 
ular bosons and excess fermions such that the molecu- 
lar bosons can dissociate due to boson-boson or boson- 
fermion interactions, but be in chemical equilibrium with 
excess fermions. When dissociation of molecular bosons 
is included, the system is no longer a binary mixture of 
molecular bosons and excess fermions, but a ternary mix- 
ture of molecular bosons, dissociated bosons (tl^T + I); 
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and excess fermions, or effectively a ternary mixture of 
molecular bosons, and f- and j-type fermions. In the 
case of ternary mixtures, there can be a large number 
of phase separated regimes. If we confine our discus- 
sion to the equillibrium of a maximum of two phases of 
this ternary mixture then several other situations can 
be encountered. For example, when j-type fermions are 
in excess, a possible sequence of phases for fixed popu- 
lation imbalance P and increasing scattering parameter 
1/{Kf<if) i s: (1) Normal phase (N) of partially polarized 
fermions — > (2) mixture of molecular bosons and T _ type 
fermions phase separated from partially polarized normal 
fermions — > (3) molecular bosons phase separated from 
excess j-type fermions — > (4) mixture of molecular bosons 
and |-type excess fermions phase separated from |-type 
excess fermions — > (5) coexistence of molecular bosons 
and |-type excess fermions. 

Therefore, as long as Fermi-Fermi mixtures can be re- 
garded as a binary mixture of non-dissociated molecu- 
lar bosons and excess fermions, the expressions given 
in Eqs. (|24|) and (j25|) may be used as a guide for the 
boundaries between phase separated (non-uniform) and 
the mixed (uniform) phases for any mixture of fermions. 
In particular, Eqs. (|24|) and ([25]) may serve as estimators 
for the phase boundaries of equal or unequal mass Fermi- 
Fermi mixtures with population imbalance, as discussed 
next. 



C. Phase Diagrams of Fermi-Fermi Mixtures in the 
Strong Attraction Limit 

Among many possibilities of Fermi- Fermi mixtures (see 
Table |T]), we focus our analysis on population imbalanced 
mixtures of 6 Li or 40 K atoms where = mi, and 6 Li 
and 40 K atoms where m-\ = 0.15™^. 

In Fig. [3l we show phase diagrams of population im- 
balance P and scattering parameter 1/ '(Kfclf) for equal 
mass mixtures when rn| = mi, and for unequal mass 
mixtures when mj = 0.15m|. In these diagrams, we 
choose | (I) to label lighter (heavier) fermions such that 
lighter (heavier) fermions are in excess when P > 
(P < 0). Although these diagrams are strictly valid in 
the strong attraction limit when 1/(Kfcif) ^ 1, they 
are qualitatively correct when 1/(Kfcif) > 1 or as long 
as the molecular bosons are not dissociated. In the 
later case, the system may be approximately described 
as a ternary mixture of molecular bosons, f- and |-type 
fermions and many other phases are possible, as discussed 
in Sec. HVBl 

In these figures, we show the following phases: (i) the 
normal (N) phase corresponding to balanced (P = 0) or 
imbalanced (P ^ 0) mixture of unpaired |- or j-type 
fermions; (ii) uniform superfluid (U) phase where paired 
(molecular bosons) and unpaired fermions coexist; and 
(iii) phase separated (PS) non-uniform superfluid phases. 
The PS(1) region labels phase separation between pure 
excess fermions and superfluid molecular bosons, while 
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FIG. 3: Phase diagram of population imbalance P — {n-\ — 
n l)/( n T + n i) vers us scattering parameter 1/(Kfclf) for (a) 
equal masses when m-f = mj., and (b) unequal masses when 
77i| = 0.15m,| in the BEC side. We show normal (N), uniform 
superfluid (U), and phase separated non-uniform superfluid 
phases PS (1) and PS (2). 

the PS (2) region labels phase separation between pure 
excess fermions, and a mixture of excess fermions and 
superfluid molecular bosons. The phase boundary be- 
tween U and PS(2) phases is determined from Eq. (|24|) . 
and the phase boundary between PS(2) and PS(1) phases 
is determined from Eq. (|25p . For a fixed mass anisotropy, 
when \P\ is large, we find phase transitions from PS(1) 
— > PS(2) — > U as the interaction strength 1/(Kfo.f) 
increases. However, when \P\ is very small, we find a 
phase transition directly from the PS(1) to the U phase 
as 1/(Kfclf) increases. 

We would like to remark in passing that the phase dia- 
grams for mixtures of 6 Li and 87 Sr or 6 Li and 171 Yb with 
mass ratios of rriLi/msr ~ 0.068 and mLi/rayb ~ 0.035, 
which are below the critical ratio m-f /m^ w 0.073 for the 
emergence of Efimov (three body bound) states are much 
richer, since phase separation and coexistence phases in- 
volving Efimov states (trimers), molecular bosons and 
excess fermions are also present. 

It is also important to emphasize that since we use the 
exact boson-boson and boson-fermion scattering lengths, 
our phase diagrams in the strong attraction limit al- 
ready include fluctuation corrections beyond the Born 
approximation. For comparison, the corresponding phase 
diagrams within the Born approximation are described 
in Fig. 01 where the phase boundaries in the popula- 
tion imbalance P versus scattering parameter l/^Kpap) 
plane are shown for equal (mj = and unequal 
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(rn| = 0.15mjJ mass mixtures. A direct comparison of 
Figs. [3] and S] shows that the results beyond the Born ap- 
proximation a re q uantitatively different from the saddle- 
point results [U H3| in the equal mass case. These 
quantitative differences become significantly large for un- 
equal mass mixtures when heavier fermions are in ex- 
cess [H, H3, IH, [sij due to the large sensitivity of the 
exact scattering parameters on the mass ratio m j /to ^ as 
shown in Fig. [2] and Table |TJ However, the same phases 
are present in both cases, indicating that the Born ap- 
proximation captures the basic qualitative features, but 
fails to produce the phase boundaries quantitatively. 

Lastly, we would like to point out the presence of sev- 
eral triple points in the phase diagrams shown in Figs. [3] 
and [31 Along the \P\ = 1 lines, we find several triple 
points as l/(K F a F ) increases where the fully polarized 
normal phases (P = ±f ) merge with (i) the partially po- 
larized normal (N) and the PS(1) phase; or with (ii) the 
PS(1) and PS(2) phases; or with (iii) the PS(2) and U 
phases. Furthermore, there is also an additional triple 
point that occurs for small \P\ (iv) where the phases 
PS(f), PS(2) and U meet. The precise locations of these 
triple points can be obtained for any mass ratio and scat- 
tering parameter from Eqs. and (|2l)|) using the equal 
sign (=) condition. The triple point for case (ii) can be 
obtained by setting \P\ = 1 in Eq. (|23|) . and the triple 
point for case (iii) can be obtained by setting \P\ = 1 in 
Eq. (|24[) . Finally, the triple point for case (iv) can be ob- 
tained by using the equal sign (—) condition of Eqs. (|24|) 
and (f2"5|) and by solving the two equations simultane- 
ously. 

Having analyzed the phase diagrams for non-trapped 
continuous systems, we discuss next the effects of the 
trapping potential which are necessary to understand ex- 
periments involving ultracold Fermi- Fermi mixtures. 



V. TRAPPED BOSE-FERMI MIXTURES AT 
ZERO AND FINITE TEMPERATURES 

In this section, we use again the simpler description 
of the effective Bose-Fermi mixture to describe trapped 
Fermi-Fermi mixtures in the strong attraction limit. For 
this purpose, we present first the theory of trapped Bose- 
Fermi mixtures at zero and finite temperatures, and then 
discuss the density profiles of trapped Fermi-Fermi mix- 
tures in the strong attraction limit using the relation be- 
tween the two systems described in Sec. IIV1 



A. Bose-Fermi mixtures 

The Hamiltonian density for a Bose-Fermi mixture in 
an external potential can be written as 

H(r) = K F (r) + V F {r)n F (r) + K B (r) + V B (r)n B (r) 
+ U BF n B (r)n F {r) + (U BB /2)n B (r)n B (r), (26) 
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FIG. 4: Phase diagram of population imbalance P — (n| — 
n l)/( n T + n l) versus scattering parameter 1/(Kfclf) in the 
Born approximation for (a) equal masses when m-f = mj., and 
(b) unequal masses when m-f = 0.15m.j in the BEC side. We 
show normal (N), uniform superfluid (U), and phase separated 
non-uniform superfluid phases PS(f) and PS(2). 

where Ki(r) = ipj (r){—V 2 / (2rrii) — ^^(r) represents 
the kinetic and chemical potential terms for fermions 
(i = F) or bosons (i = B) in a single hyperfine state. 
Here, rij(r) = ipj (r)^(r) represent the local density op- 
erators, and U BF and U BB represent the boson- fermion 
and boson-boson interaction. The single-hyperfine-state 
fermions are non-interacting, but obey the Pauli exclu- 
sion principle. 

For simplicity, we approximate the trapping potential 
by an isotropic harmonic function where the potential 
energy is 

VHr) = \a ir 2 . (27) 

Here, a* = rrnw 2 is proportional to the trapping fre- 
quency of bosons [i = B) or fermions (i = F), which 
is typically different for each kind of atom. Since the 
potential and the interactions are isotropic, the effective 
potentials and densities depend only on r = \r\. 

In the presence of such trapping potentials, the bosons 
and fermions feel the effective potentials 

Vs.effW = V B (r)+2U BB n B (r) + U BF n F (r), (28) 
VF.efi(r) - V F (r) + U BF n B (r), (29) 

respectively, where n F (r) is the local density of fermions, 
and 

n B (r) = n c (r) + n NC {r). (30) 
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is the total local density of bosons. Here, nc{r) and 
nNc{ r ) are the density of condensed and non-condensed 
bosons, respectively. 

The number of condensed bosons is determined from 
the Gross-Pitaevskii equation leading to 



n c (r) 



- Vb{t) - 2U B Bn N c(r) - Ubf^f ( r ) 
Ubb 



_ (31) 

within the Thomas- Fermi approximation (TFA), where 
the kinetic energy of the bosons is neglected. This 
relation is valid when the condition /is — V B {r) — 
2UBBn-Nc( r ) ~ Ubf^f^) > is satisfied, otherwise, 
n c (r) = 0. 

In our analysis, we treat both the non-condensed 
bosons and the fermions as ideal gases [4J] in effective 
potentials given in Eqs. and (f2"9")l . In addition, we 
assume that non-condensed bosons are in thermal equi- 
librium with condensed bosons at the same chemical po- 
tential fiB- Within these approximations, the density of 
non-condensed bosons and fermions are given by 



1 



n N c(r) = ^ & [ £ k,-B - Ms + V B ^s{r)\, (32) 

k 



n F ( 



W = ^E/[ £ V-^ + ^9WL (33) 



where b(x) = l/[cxp(a;/T) — 1] is the Bose, and f(x) — 
l/[exp(x/T) + 1] is the Fermi distribution. Here, ek,i = 



k| /(2nii) is the kinetic energy of bosons 



B) or 



fermions (i = F). 

Notice that, at zero temperature, all bosons condense 
such that ?i jvc ( r ) — 0, and riB{r) = nc{r), leading to 



n B {r) 



VB - VM/) - Ubfit-f (x) 
Ubb 



(34) 



, x {2m F \pL F - V F (r) - U B Fn B (r)}} 3/2 
n F (r) = — 2 ,(35) 

for the densities of bosons and fermions, respectively [3pj ] . 
The first expression is valid when the condition /is — 
Vs{r) — UsFTiFir) > is satisfied, otherwise, riB{r) = 0. 
The second expression is valid when the condition fi F — 
Vp{r) — UsFnB(r) > is satisfied, otherwise, n F {r) = 0. 

The chemical potentials of bosons and fermions are 
determined from fixing the total number of bosons and 
the number of fermions, independently, as follows 



N F = 



d 3 rn B (r), 
d 3 rn F {r), 



(36) 
(37) 



where the integration is over all space. Therefore, in 
order to find the density profiles for condensed and non- 
condensed bosons, as wells as for fermions, we need to 
solve Eqs. ([36]) and (|37| for hb and /i F self-consistently. 
Next, we discuss the density profiles of Fermi- Fermi mix- 
tures in the strong attraction limit using the effective 
Bose-Fermi description presented. 



B. Fermi- Fermi Mixtures 
in the Strong Attraction Limit 

To make the connection between Bose-Fermi mixtures 
and population imbalanced Fermi-Fermi mixtures in the 
strong attraction limit, we identify F = {| or |} as the 
excess fermions. This identification leads to the density 
of excess fermions (tie) and molecular bosons {n B ) given 
by 



n E (r) = n F (r) - n- F (r) = |n T (r) - nj.(r)|, (38) 

(39) 



, , n(r) - n E (r) 
nB{r) = =n_ F (r), 



respectively, where n(r) — n^{r) + n±(r) is the total den- 
sity of f- and j-type fermions. Here, we use (— f ) =J. and 
vice versa. For instance, if F =f fermions are in excess, 
the density of excess fermions and molecular bosons are 
nE{f) = Uf(r) — ni(r) and ns(r) = [n(r) — n E (r)]/2 = 
ni(r), respectively, such that all |-type fermions are 
paired with some of the f-type fermions to form molec- 
ular bosons, but there are T-type fermions left unpaired. 
In addition, we identify as — a^+a^, where a a = m cr w' 1 C! 
is different for different atoms. 

As an example, we look at the equal mass case m F = 
TOf = mi = m and m/j = m-f + = 2m, and solve 
the self-consistent relations Eqs. (|3€> j) and ([57)) for a two- 
hyperfine-state mixture of 6 Li only or K only. In our 
numerical analysis, we choose a F = = = a and 
olb = 2a, scattering parameter l/(K F a F ) — 2 and popu- 
lation imbalance parameter P = 0.5 such that n-\ = 3n^. 
We define a characteristic energy E F = K F /(2m) in 
terms of the wavevector K F and the fermion mass m 
and assume it to be the known Fermi energy of a gas of 
non interacting fermions in the trapping potential Vp (r) . 
We scale the radial distance with the effective Thomas- 
Fermi (TF) radius R TF defined by E F = aR\ F j2. With 
this identification, the total number of trapped fermions 
N = (E F /w F ) 3 /3, with wp — y/a/m, can be rewritten 
as N — KpR 3 r, F /24 in terms of K F and Rtf- 

Since we are interested also in the temperature depen- 
dence of the density profiles, we recall that the critical 
temperature for Bose-Einstein condensation of a non- 
interacting harmonically trapped Bose gas is Tbec — 
wb[Nb I 'C(3)] 1 ^ a ', where wb = \J 1 olbIvhb- In our Fermi- 
Fermi mixture the number of bosons is Nb = N(l — 
\P\)/2 expressed in terms of the total number of fermions 



N 



N t 



N± and population imbalance P = (N-\ 



N±)/N. Using the expression for Nb, the expression of 
N in terms of Ep and w F , and that wb — w f for equal 
masses, we find 



T. 



BBC — 



6C(3) 



E F ^0.518(1 -\P\) 1/3 E F , (40) 



for the critical BEC temperature, which is valid when 
\/{Kpap) — > oo. Here, (,{x) is the Zeta function and 
C(3) w 1.202. Therefore, for P = 0.5, we obtain T BE c ~ 
OAlEp. Notice that Tbec for population imbalance P 
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FIG. 5: Density (in units of K F ) of condensed (nc) and non- 
condensed (njvc) molecular bosons, and excess (tie) fermions 
versus trap radius t/Rtf is shown for (a) T = 0, (b) T — 
0.2E F , (c) T = 0.35E F , and (d) T = 0A1E F . Here, m T = 
my, the population imbalance parameter is P — 0.5 and the 
scattering length parameter is l/(K F a F ) = 2. 



given in Eq. P0")) is a generalization of the results for 
equal populations [H, [46[ • 

In Fig. [5l we show the density (in units of Kp) of con- 
densed (nc) and non-condensed (njsrc) molecular bosons, 
and excess (tie) fermions as a function trap radius r / Rtf 
for four temperatures: (a) T = 0, (b) T = 0.2£ F , (c) 



T = 0.35E F , and (d) T = 0AIE F . At zero tempera- 
ture (T = 0), as shown in Fig. ODJa), we find that all of 
the molecular bosons are condensed, and that they are 
concentrated close to the center of the trap. In contrast, 
the majority of excess fermions are pushed away from the 
center towards the edges of the trap due to the repulsive 
boson-fermion interaction and the high concentration of 
condensed molecular bosons. Therefore, there is a clear 
indication of phase separation between molecular bosons 
and excess fermions. When the temperature is increased 
to T = 0.2E F shown in Fig. O^b), some of the molecular 
bosons are not condensed. These non-condensed molec- 
ular bosons are also pushed away from the center to- 
wards the edges of the trap just like the excess fermions. 
Further increase in temperature increases (decreases) the 
number of non-condensed (condensed) molecular bosons 
as can be seen in Fig. 03c). For temperatures close to 
Tbec and above, all of the molecular bosons become non- 
condensed as shown in Fig. [5Jd) having a Gaussian-like 
density distribution. Similarly, the excess fermions also 
have Gaussian-like density distribution for temperatures 
at Tbec and above due to the absence of the condensate. 

Therefore, at zero temperature, we find that the trap- 
ping potential tends to favor phase separation into a 
PS(l)-rich phase where regions of almost pure fermions 
and almost pure bosons are separated. However, at finite 
temperatures, the system develops a PS(2)-rich phase 
where regions of almost pure fermions and of almost fully 
mixed bosons and fermions are separated. The region of 
coexistence of bosons and fermions can be further bro- 
ken down into a domain of coexisting excess fermions 
with condensed and non-condensed bosons, and into a 
domain of coexisting excess fermions and non-condensed 
bosons as can be seen in Fig.[5jd). Again, if the molecular 
bosons are allowed to dissociate, the system can be de- 
scribed by a ternary mixture as discussed in Sec. IIV Bl and 
the phase diagram can be even richer, especially closer 
to the unitarity. 

In Fig. [HI we show the density (in units of Kp) of f- 
(rif) and |-type (njj fermions as a function of radius 
r/R TF for four temperatures: (a) T = 0, (b) T = 0.2E F , 
(c) T = 0.35E F , and (d) T = 0A1E F . At zero tem- 
perature, as shown in Fig.^a), we find that the density 
of T- and |-type fermions are similar close to the cen- 
ter of the trapping potential, while some of the excess- 
type fermions are close to the edges. When the temper- 
ature increases to T = 0.2E F shown in Fig. Hfb) or to 
T = 0.35E F shown in Fig.[6{c), the density of \- and J.- 
type fermions become different at the center of the trap. 
In addition, both f- and J.-type fermions exist towards 
the edges. At temperatures close to Tbec and above, the 
density profiles of f- and j-type fermions have the stan- 
dard shapes of weakly interacting trapped Fermi gases. 

In this section, we have shown that the effective Bose- 
Fermi description of Fermi-Fermi mixtures is applica- 
ble in the strong attraction limit, and thus can pro- 
vide good quantitative comparisons to experiments in the 
same regime, since the exact boson-fermion and boson- 
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FIG. 6: Density (in units of Kp) of j- (n-f) and |-type (n^) 
fermions versus trap radius v/Rtf is shown for (a) T = 0, 
(b) T = 0.2E F , (c) T = Q.35E F , and (d) T = QAIEf. Here, 
m-f = mj., population imbalance parameter is P — 0.5 and 
scattering length parameter is 1/{Kfclf) = 2. 



boson scattering parameters were used to obtain the 
phase diagrams and density profiles. Thus, we think 
that it is particularly important to perform experiments 
for different population imbalances in the strong attrac- 
tion limit, where the theory is simple. The situation is 
somewhat more complicated near unitarity where quan- 
titative comparisons between theory and experiment are 



more difficult. Furthermore, there are also some differ- 
ences between the experimental studies of MIT 0, Q and 
Rice [1, [l(| groups performed near unitarity, since the 
shapes of their traps and the number of trapped atoms 
are quite different. 

Having concluded the analysis of the effects of trapping 
potentials on Fermi-Fermi mixtures in the strong attrac- 
tion limit, next we give a summary of our conclusions. 



VI. CONCLUSIONS 

In summary, we used the effective Bose-Fermi mixture 
description to obtain the phase diagrams of Fermi- Fermi 
mixtures with equal or unequal masses and equal or un- 
equal populations in the strong attraction limit. For 
this purpose, we analyzed first the exact boson-fermion 
and boson-boson scattering lengths as a function of mass 
anisotropy, and then we constructed the phase diagrams 
of Fermi-Fermi mixtures in the BEC regime. 

We showed that three-dimensional non-trapped 
fermion mixtures with population imbalance exhibit 
phase separation in addition to the normal polarized mix- 
ture of fermions and uniform mixture of superfluid and 
excess fermions. In the BEC regime, we found two dif- 
ferent non-uniform phase separated states: PS(1), where 
there is phase separation between pure unpaired (ex- 
cess) and pure paired fermions (molecular bosons); and 
PS (2), where there is phase separation between pure ex- 
cess fermions and a mixture of excess fermions and molec- 
ular bosons. For equal mass mixtures, our results for the 
phase boundaries are quantitatively different from pre- 
vious saddle-point results, and these quantitative differ- 
ences become more pronounced for unequal mass mix- 
tures when heavier fermions are in excess indicating the 
importance of taking into account scattering processes 
beyond the Born approximation. 

We also discussed the effects of trapping potentials 
on the density profiles of condensed and non-condensed 
molecular bosons, and excess fermions at zero and finite 
temperatures. At zero temperature, we found that al- 
most all of the condensed bosons are at the center of the 
trap, while the excess fermions are pushed to the edges 
due to the repulsive boson-fermion interactions. At finite 
temperatures, we found that non-condensed pairs and 
excess fermions are created at the center of the trap at 
the expense of an overall reduction of condensed bosons. 
Finally, at temperatures above the BEC temperature, 
the number of condensed bosons vanish, and the system 
becomes a mixture of weakly interacting non-condensed 
bosons and excess fermions. Finally, we discussed that 
our findings can provide good quantitative comparisons 
to experiments performed in the same regime of validity 
of the theory (BEC regime), since the boson-fermion and 
boson-boson scattering parameters that enter our calcu- 
lations are exact in the dilute limit. 

Lastly, we think that it is important to perform experi- 
ments with Fermi-Fermi mixtures in the strong attraction 
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limit (BEC regime) where the theoretical description is observed, and directly compared with the theory, 
simple. In this limit, additional superfluid and normal 

phases and richer density profiles proposed here can be We thank NSF (DMR-0304380) for support. 
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